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INTRODUCTION 


An  investigation  is  being  conducted  at  the  Harry  Diamond  Laborato¬ 
ries  { HDL )  to  develop  a  Cerenkov  device  as  a  relatively  intense  source 
of  near-millimeter  radiation*  Sources  (both  pulsed  and  continuous)  that 
produce  high  power  are  desired  for  many  of  the  projected  applications  of 
near-millimeter  waves*  The  experiment  initiated  at  HDL  represents  one 
approach  to  designing  such  a  radiation  source* 

The  foundation  for  this  experiment  has  been  laid  by  Felch  et  al  at 
Dartmouth  College.1'2  Very  simply#  that  experiment  consists  of  passing 
an  intense,  relativistic  electron  beam  through  a  cylindrical  waveguide 
with  a  thin  dielectric  liner.  The  beam  current  is  pulsed,  having  a  peak 
value  of  1  to  10  A,  and  the  electrons  are  accelerated  to  relativistic 
speeds  through  the  dielectric-lined  resonator,  where  the  actual  produc¬ 
tion  of  radiation  takes  place. 

As  a  result  of  the  considerable  success  of  the  Dartmouth  device,  HDL 
has  undertaken  a  similar  experiment.  Dartmouth's  pulsed  device  has 
generated  powers  on  the  order  of  hundreds  of  watts  and  at  fairly  widely 
spaced  frequencies.1  However,  the  detailed  characterization  of  the 
output  has  proven  difficult  partially  because  of  effects  introduced  by 
the  pulsing  process.  To  overcome  these  difficulties,  a  continuous  wave 
(cw)  device  is  being  investigated  at  HDL.  The  system  will  have  a  200-kV 
electron  beam  with  a  maximum  continuous  current  of  about  250  mA.  In 
most  other  respects,  the  two  systems  will  be  similar. 

This  report  deals  with  the  design  of  the  resonating  structure.  The 
basic  principles  for  the  design  of  a  Cerenkov  power  source  have  been 
worked  out  to  a  considerable  degree  at  Dartmouth  College.2'3  Since  the 
theories  underlying  the  HDL  and  Dartmouth  experiments  are  fundamentally 
the  same,  the  Dartmouth  results  are  used  as  a  basis  for  the  design. 

Radiation  in  the  Cerenkov  device  is  generated  by  an  interaction 
between  the  electron  beam  and  the  partially  filled  waveguide.  Deriving 
the  dispersion  relation  for  a  relativistic  electron  beam  in  a  partially 
filled  waveguide  reveals  that  the  phase  velocities  of  modes  that  can 
propagate  in  the  guide  are  all  superluminal,  that  is,  greater  than  the 


1/C.  L  .  Felch ,  K*  O.  Busby,  R .  W .  Layman ,  and  J.  E  .  Walsh ,  Cerenkov 
and  Cerenkov- Raman  Masers:  Experiments ,  Physics  of  Quantum  Electronics , 

7_  (1980),  2  55-300  and  301-322. 

2/C •  L.  Felch,  Cerenkov  Radiation  in  Dielectric-Loaded  Waveguides , 

Doctoral  Dissertation,  Dartmouth  College,  Dartmouth ,  NH  (1980). 

3M.  Johnson,  Cerenkov  Interaction  in  a  Waveguide  Partially  Filled 
with  Dielectric ,  Master's  Thesis,  Dartmouth  College,  Dartmouth,  NH 
(1978). 

\ 
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speed  of  light  in  the  dielectric.  By  using  a  beam  with  a  particle 
velocity  that  also  is  superluminal ,  it  is  possible  to  couple  the  beam  to 
a  given  waveguide  mode,  thus  driving  an  instability  with  a  sizable 
growth  rate. 

The  problem  discussed  here  is  the  use  of  these  theoretical  consider¬ 
ations,  approximate  though  they  are,  in  the  design  of  a  resonating 
structure  for  the  proposed  power  source.  The  basic  parameters  are  beam 
current,  beam  velocity  (voltage),  dielectric  constant  of  the  liner, 
radius  of  the  waveguide,  and  thickness  of  the  liner.  Since  this  is 
intended  as  a  near-millimeter  wave  source,  an  output  frequency  of  120 
GHz  will  be  assumed. 

The  remainder  of  this  report  begins  with  a  brief  section  on  the 
theoretical  foundation.  This  is  followed  by  an  analysis  that  first  uses 
the  no-beam  dispersion  relation  derived  by  Felch2  to  calculate  modes  for 
parameters  relevant  to  our  experiment.  Next,  our  analysis  continues  by 
using  a  program  that  we  have  written  to  solve  for  higher  order  modes  in 
the  resonator  to  consider  the  possibility  of  coupling  the  beam  to  higher 
modes.  Finally,  there  is  a  summary  of  what  can  be  concluded  from  the 
analysis • 


2 .  THEORY 

The  derivation  outlined  here  parallels  the  work  of  Felch.2  The 
dispersion  relation  for  this  system  is  obtained  by  using  a  linear  per¬ 
turbation  method  and  by  making  the  following  assumptions :  ( 1 )  the 

dielectric-lined  waveguide  is  infinitely  long,  (2)  the  metal  portions  of 
the  guide  are  perfectly  conducting,  and  (3)  the  magnetic  waveguide  field 
is  sufficiently  strong  to  make  the  transverse  velocity  of  the  electrons 
in  our  beam  negligible.  The  frequency  and  the  magnetic  field  proposed 
for  this  experiment  justify  assumptions  (2)  and  (3). 

The  geometry  being  considered  is  shown  in  figure  1,  in  which  the 
center  region  is  evacuated  and  entirely  filled  by  the  electron  beam. 


z 


Figure  1.  Geometry  of  beam-driven 
cylindrical  waveguide  with 
dielectric  liner. 


2fC.  L.  Felch,  Cerenkov  Radiation  in  Dielec  trie -Loaded  Waveguides , 
Doctoral  Dissertation,  Dartmouth  College,  Dartmouth,  NH  (1980). 
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The  following  form  of  Maxwell's  equations  is  used: 


•  D  =  4  Tip  =  -47ine  , 


?  •  B  *  0  , 


±  t  4tt  t  1  3D 
V  x  B  -  —  3  +  -  —  r 

c  c  3t 


±  £  1  3B 

V  x  E  =  -  —  — 
c  3t 


Also ,  the  Lorentz  force  and  continuity  equations  are  used 


md(yv)  +  e  ■*  x 

— -J =  -eE  -  —  v  x  B  , 

dt  c 


+  $  •  (nv)  =  0  , 


where 

*  z  gE,  a  <  r  <  b, 

D  ={> 

E,  0  <  r  <  a , 

g  =  dielectric  constant  of  liner, 

E  =  electric  field  strength, 

n  =  electron  density  in  beam, 

-e  =  electronic  charge, 

-► 

B  =  magnetic  induction, 

c  =  velocity  of  light  in  vacuum, 

>  + 
j  **  current  density  =  -env, 

v  =  velocity  of  electrons  in  beam, 

m  =  electronic  mass, 

Y  =  1/(1  -  v2/c*)l/2. 
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In  the  general  case  of  a  partially  filled  waveguide,  transverse 
magnetic  (TM)  ,  transverse  electric  (TE) ,  and  hybrid  modes  can  propagate 
in  the  waveguide.4'5  The  discussion  is  restricted  to  TM  modes  since  a 
complete  analysis  reveals  that  the  beam  will  couple  to  the  z-component 
of  the  electric  field.  Further,  because  of  the  cylindrical  symmetry  of 
the  guide,  only  TMgn  modes  need  be  considered. 

The  field  components  of  the  TM^n  mode  are 


E  =  E  (r)eikz“ia,te  +  E  (r)eikz“iwte 
r  r  z  z 


(7) 


where 


,ikz-i(jdtf 

e8 


+  Be 
0  z 


(8) 


k  =  guide  wave  number, 
u)  =  angular  frequency, 

Er  =  component  of  E  in  r  direction, 

Ez  =  component  of  E  in  t  direction, 

Bq  =  component  of  B  in  0  direction, 

Bq  =  external  magnetic  field. 

In  general,  these  fields  will  perturb  n  and  v  with  plane-wave 
components.  The  result  is 


n  =  nQ  +  n^rJe 


ikz-iwt 


(9) 


v  =  v  e  +  v  (r)eikz~iwte 
0  z  1  z 


(10) 


In  equations  (9)  and  (10)  and  subsequent  equations,  the  M0"  subscript  on 
n,  v,  and  j  refers  to  a  zero-order  quantity  and  the  **  1"  subscript  on  n, 
v,  and  j  refers  to  a  first-order  approximation. 

If  one  assumes  that  the  fields  and,  correspondingly,  n^(r)  and  v-j(r) 
are  small,  then  only  first-order  terms  need  be  kept  and 


4M.  Miyagi  and  S.  Nishida ,  Transmission  Characteristics  of  Dielectric 
Tube  Leaky  Waveguide ,  IEEE  Trans •  Microwave  Theory  Tech.,  MTT-28  (June 
1980),  536-541 . 

5L.  P incherle ,  Electromagnetic  Waves  in  Metal  Tubes  Filled 
Longitudinally  with  Two  Dielectrics ,  Pbys.  Rev .,  66_  (September  1944), 
118-130 . 


+  > 
j  =  nev 


=  £-riQevQ  -  -  nQev-j  (  r )  e^z-iuJtj  (11! 


jQ  -  j,(r) 


Now,  using  equations  (5)  and  (€>),  one  may  find  expressions  for  n-j(r) 
and  v^(r)/  which  when  substituted  into  equation  (11)  lead  to 


jAr)  =  - o7? - 77  Ez<r>e 

'  4tty3(u)2  -  kvo)2  2 

where  gj2  =  4ffn  e2/m  is  the  plasma  frequency  of  the  electron  beam* 

p  0 

Manipulation  of  Maxwells  equations  readily  yields 


-T^)- 


V2E  =  Hi»  t  u£e  *  , 

c2  1  c2 


the  z-component  of  which  is 


rdE  (r) 


1  A.  zvx 

r  Sr  3r 


+  qzSzEz(r)  =  0  , 


where 


2  _  go2e  t  ? 

q  ^ - k2 

^2 


=  1  - 


£Y3(o)  ~  kvQ)2 


and  equation  (12)  is  used  for  }  . 


Equation  (14)  is  Bessel’s  equation  and  has  the  general  solution 


E  (r)  =  AJ  ( qSr )  +  BY(qSr) 
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where  Jg  and  Yg  are  zero-order  Bessel  functions  of  the  first  and  second 
kinds,  respectively. 

In  the  same  way,  examination  of  the  r-component  of  equation  (13) 
yields 


3E 

lk  z 

E  =  -  — - 

r  o  3r 

q 


(  18) 


Consider  the  two  regions  in  the  waveguide  with  the  conditions 

£  =  £ 


and 


w2  =  0 
P 

£  =  1 
U)2  *  0 

P 


!  . 


a  <  r  <  b 


r  <  a  • 


(  19) 


(20) 


Then 


and 


(r)  =  -  i|  [AJ^qr)  +  BYj (gr)  ]  'j 


Ez(r)  =  AJQ(qr)  +  BYQ(qr) 


a  <  r  <  b 


(21  ) 


and 


■j  If  T1  r 

E  (r)  =  -  f  A*  J  ( ipTr )  +  B’Y  (ipTr) 

r  ip  1  1  u 


E  (r)  =  A’J.(ipTr)  +  B’Y  (ipTr) 
z  u  u 


\  ,  r  <  a  ,  ( 22 ) 


where  and  Y1  are  first-order  Bessel  functions  of  the  first  and  second 
kinds,  respectively, 


T2  =  1 - £ -  ,  (23) 

Y3(u>  -  kvQ)2 

p2  i  k2  -  ,  (24) 

c2 

and  A,  B,  A',  and  B*  are  constants.  The  constants  are  evaluated  by 
applying  the  usual  boundary  conditions: 
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•  E  =  0  at  r  =  b. 

•  E  is  finite  at  r  =  0. 

•  The  tangential  component  of  E  is  continuous  at  r  =  a . 

•  The  normal  component  of  eE  is  continuous  at  r  =  a . 


Applying  these  boundary  conditions  results  in  the  following 
transcendental  relation : 


J1 (ipTa) [j0(qb)Y0(qa)  -  YQ ( qb ) JQ ( qa )  ] 


(25) 


+  |  J0(ipTa)[j()(qb)Y1(qa)  -  YQ ( qb ) J 1 ( qa )  ]  =  0 

which  is  the  dispersion  relation  from  TMQn  modes  for  this  system. 

In  the  next  section,  we  consider  this  dispersion  relation  in  the 
context  of  our  experiment. 


3.  ANALYSIS 

3.1  System  Parameters 

The  calculations  to  be  described  have  been  made  for  parameters 
similar  to  those  that  we  expect  to  have  with  our  experiment.  The  goal 
of  producing  120-GHz  (wavelength  of  2.5  mm)  radiation  encourages  the  use 
of  as  small  a  waveguide  as  is  practical.  On  the  other  hand,  beam  powers 
of  many  kilowatts  and  the  problems  of  fabrication  force  us  to  use  a 
larger  waveguide  than  would  otherwise  be  necessary  for  120-GHz  opera¬ 
tion  • 


We  deal  with  beam  velocities  no  greater  than  0.7c  (200  kV),  so  a 
relatively  large  dielectric  constant  is  desirable  to  assure  that  the 
particle  velocity  is  well  into  the  superluminal  region.  One  possibility 
is  fused  quartz  with  e  =  3.78.  The  following  are  then  a  reasonable  set 
of  initial  values  for  our  analysis  : 

•  Output  frequency  =  120  GHz. 

•  Relative  beam  velocity,  3  =  v/c  =  0.7. 

•  Inner  waveguide  radius,  a  =  0.45  cm. 

•  Dielectric  constant,  c  -  3.78. 


11 


3.2  Dispersion  Curves 

If  we  choose  a  fixed  value  for  u),  then  there  are  many  values  of 
k  that  satisfy  equation  (25).  Each  one  of  these  roots  corresponds  to 
the  wave  number  of  a  mode  that  will  propagate  in  the  waveguide  at  fre¬ 
quency  0).  The  values  of  k  corresponding  to  the  fundamental  mode,  TMq  ^ , 
for  varying  u)  have  been  found  by  Johnson.3  However,  solving  equation 
(25)  in  this  way  is  fairly  difficult  since  k  is  generally  complex. 
Fortunately,  as  demonstrated  by  Johnson,  even  for  beam  currents  much 
higher  than  those  that  we  intend  to  produce,  the  no-beam  approximation 
is  very  close  to  the  exact  solution.  This  is  the  approximation  used  by 
Felch2  as  well. 

The  no-beam  approximation  results  from  setting  w  =  0  in  equa¬ 
tion  (25).  Thus,  it  reduces  to  ^ 


P 


J 1  (  ipa  )  [  JQ  ( qb  )  YQ  ( qa  ) 


Yo(qb)Jo(qa) J 


(26) 


+  ^  JQ(ipa) [jQ(qb)Y1 (qa)  -  YQ ( qb) J ^ ( qa )  ]  =  0 

A  convenient  way  of  examining  this  relation  is  to  plot  w/c  as  a 
function  of  k  while  varying  the  different  parameters.  In  figure  2, 
these  are  plotted  for  different  dielectric  thicknesses  for  the  TMq  ^ 
mode.  Since  the  beam-to-dielectric  distance  is  critical,  we  have  fixed 
the  inner  radius  at  a  *  0.45  cm  and  varied  b.  For  b/a  =  1.0,  the  dis¬ 
persion  curve  asymptotically  approaches  that  for  a  plane  wave  in  a 
vacuum.  This  is  the  expected  result  for  any  empty  waveguide.  As  we 
consider  dielectric  liners  of  increasing  thicknesses,  the  phase  velocity 
of  the  modes  approaches  the  speed  of  light  in  the  dielectric. 

The  curve  w/ck  =  3  “  0.7  represents  the  "beam  curve,"  that  is, 
the  dispersion  curve  of  the  electron  beam.  We  note  where  the  different 
waveguide  curves  cross  the  beam  curve.  At  this  frequency,  called  the 
synchronous  frequency,  the  maximum  transfer  of  energy  from  the  beam  to 
the  waveguide  mode  should  occur.3  Clearly,  as  the  dielectric  gets 
thinner,  the  synchronous  frequency  becomes  higher.  It  is  easy  to  see 
from  figure  2  that  we  may  vary  the  synchronous  frequency  also  by  varying 
the  beam  velocity. 

.  L.  Felch,  Cerenkov  Radiation  in  Dielectric-Loaded  Waveguides, 
Doctoral  Dissertation ,  Dartmouth  College ,  Dartmouth ,  NH  (1980)* 

3M.  Johnson,  Cerenkov  Interaction  in  a  Waveguide  Partially  Filled 
with  Dielectric ,  Waster's  Thesis,  Dartmouth  College,  Dartmouth ,  NH 
(1978). 
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HAM  CURVE 


Figure  2.  Behavior  of  waveguide 
dispersion  relation  for  different 
dielectric  liner  thicknesses,  TM01 
mode:  3  =  0*7,  e  =  3.78,  and 

a  =  0.45  cm. 


Figure  3  shows  the  case  where  the  beam  couples  to  the  waveguide 
at  120  GHz.  This  coupling  requires  a  thickness  of  b/a  -  1.085,  which 
corresponds  to  a  wall  thickness  of  0.3825  mm  (0.015  in.). 


Figure  3.  Behavior  of  waveguide 
dispersion  relation  for  b/a  =  1.085, 
TMQ1  mode:  6  =  0-7,  e  =  3.78, 
and  a  *  0.45  cm. 


These  curves  were  all  generated  by  using  e  =  3.78?  however,  as 
shown  in  figure  4,  we  may  use  the  dielectric  constant  to  choose  a  cou¬ 
pling  frequency  if  various  dielectrics  are  available  to  us.  Comparison 
of  figures  2  and  4  indicates  that  the  dispersion  curves  are  more  sensi¬ 
tive  to  the  dielectric  thickness  than  to  the  dielectic  constant. 


Figure  4* 
dispersion 
dielectric 
3  =  0.7,  a 


Behavior  of  waveguide 
relation  for  different 
constants,  TMq i  mode: 

=  0.45  cm,  and  b/a  =  1. 


1  . 


A  fourth  way  to  change  the  characteristics  of  the  waveguide  is 
to  vary  the  inner  radius,  a,  of  the  guide.  Figure  5  illustrates  this 
variance  for  several  values  of  a  from  0.35  to  0.55  cm.  Again,  the 
curves  are  less  sensitive  to  the  inner  radius  than  to  the  liner  thick¬ 
ness.  Clearly,  there  are  numerous  ways  to  select  a  preferred  frequency 
for  coupling.  To  a  large  extent,  the  method  we  use  depends  on  the 
detailed  mechanical  design  of  the  system. 


Figure  5.  Behavior  of  waveguide 
dispersion  relation  for  different 
inner  radii,  TMq ^  mode:  3  =  0.7, 
e  =  3.78,  and  b/a  =  1.1. 


3 • 3  Higher  Order  Modes 

The  program  that  generated  figures  2  to  5  was  an  adaptation  of 
the  one  written  by  Felch2  and  calculates  only  the  fundamental  (TMq-j) 
mode.  It  is  conceivable  that  we  may  wish  to  couple  to  a  higher  order 

•  L •  Felch,  Cerenkov  Radiation  in  Dielectric-Loaded  Waveguides , 
Doctoral  Dissertation ,  Dartmouth  College ,  Dartmouth,  NH  (1980). 
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mode.  Fbr  example ,  the  situation  represented  in  figure  3  requires  a 
dielectric  wall  of  0.3825  mm,  whereas  coupling  to  a  higher  mode  would 
allow  a  larger  thickness,  which  would  be  easier  to  fabricate.  We  are 
interested  also  in  the  separation  of  the  different  modes;  therefore,  it 
was  instructive  to  write  a  program  to  calculate  the  dispersion  relations 
for  the  higher  order  TMgn  modes.  In  this  calculation,  we  assume  that 
the  previous  no-beam  approximation  remains  valid  for  the  higher  modes. 

Figure  6  shows  the  higher  modes  for  b/a  =  1.085  and  e  =  3.78. 
As  in  figure  3,  the  synchronous  frequency  for  TMq  ^  is  about  120  GHz  ( co/c 
=  25.13).  In  figure  7,  b/a  =  1.3  and  the  TMQ2  mode  couples  to  the  beam 
at  the  desired  frequency.  Similarly,  figure  8  shows  this  condition  for 
TMq^  with  b/a  =  1.5.  Thus,  we  can  couple  to  higher  modes  by  using 
thicker  liners.  Choosing  higher  modes  is  feasible  also  by  varying  the 
beam  velocity  or  the  dielectric  constant  as  we  have  demonstrated  for  the 
fundamental  mode. 


Figure  6.  Dispersion  relations  for 
higher  order  TMQn  modes:  0  =  0.7, 
e  =  3.78,  a  =  0.45  cm,  and  b/a  = 

1.085. 


Figure  7.  Dispersion  relations  for 
higher  order  TMpn  modes:  0  =  0.7, 
e  =  3.78,  a  =  0.45  cm,  and  b/a  = 

1.3. 


When  coupling  to  higher  modes,  one  has  to  consider  the  possibil¬ 
ity  of  mode  competition  from  lower  modes  as  well  as  from  modes  that  are 
higher  still.  One  can  see  from  figure  7  that,  while  coupling  to  TMQ2  at 
120  GHz,  we  couple  also  to  TMq^  at  38  GHz.  Similarly,  figure  8  demon¬ 
strates  that  coupling  to  TMg^  at  our  desired  frequency  results  in  cou¬ 
pling  also  to  TMQ1  at  23  GHz  and  to  TMQ2  at  71  GHz.  Although  the  fig¬ 
ures  presented  here  do  not  illustrate  coupling  to  higher  modes  (TMq^, 
TMqs /  etc.),  this  coupling  can  occur  and  should  be  considered.  There¬ 
fore,  one  needs  to  examine  the  competition  between  these  different  modes 
to  determine  the  practicality  of  coupling  to  a  higher  order  mode. 
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Figure  8.  Dispersion  relations  for 
higher  order  TMQn  modes:  3  =  0.7, 
e  -  3.78,  a  =  0.45  cm,  and 
b/a  =  1.5. 


4.  CONCLUSIONS 

Presented  here  is  an  approach  to  designing  the  resonator  for  the 
Cerenkov  device  based  on  the  theoretical  dispersion  relation  and  some 
physical  restrictions  of  the  system.  It  has  been  demonstrated  that  by 
using  the  proper  dimensions,  dielectric,  and  beam  velocity,  we  can 
reasonably  expect  to  couple  the  beam  to  low  order  TM  modes  at  the  de¬ 
sired  frequency. 

q 

When  solving  the  exact  dispersion  relation,  Johnson0  suggests  that 
the  imaginary  roots  contain  information  about  the  growth  rates  of  the 
instability  being  driven  by  the  beam.  Fortunately,  his  results  indicate 
that  the  parameters  given  in  section  3.1  should  produce  relatively  high 
growth  rates,  provided  that  the  beam  energy  is  sufficient.  Further 
analysis  of  the  growth  rates  should  yield  useful  information  regarding 
mode  competition  when  coupling  to  higher  order  modes  as  well  as  helping 
to  maximize  the  coupling  for  the  fundamental  mode. 

This  model  of  our  resonating  structure  is  not  exact.  Numerous 
approximations  were  made,  including  a  very  narrow  velocity  distribution 
in  the  electron  beam  and  an  infinitely  strong  and  straight  magnetic 
field  for  the  beam.  But  perhaps  the  most  unrealistic  approximation  was 
the  assumption  of  an  infinitely  long  resonator.  For  these  calculations 
to  be  representative  of  the  experimental  arrangement,  one  would  probably 
need  to  do  some  impedance  matching  at  both  ends  of  the  dielectric  liner 
to  minimize  reflection  of  the  radiation  at  these  areas.  Nevertheless, 
this  analysis  suggests  the  parameters  to  use  in  the  initial  fabrication 
of  the  resonator. 


%.  Johnson ,  Cerenkov  Interaction  in  a  Waveguide  Partially  Filled 
with  Dielectric ,  Master9 s  Thesis ,  Dartmouth  College ,  Dartmouth,  NH 
(1978). 
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